A "HORIZONTAL" HYPER DIFFUSION 3-D THERMOCLINE PLANETARY 
GEOSTROPHIC MODEL: WELL-POSEDNESS AND LONG TIME BEHAVIOR 
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Abstract. In this paper we study a three dimensional thermocline planetary geostrophic "horizontal" 
hyper— diffusion model of the gyre-scale midlatitude ocean. We show the global existence and uniqueness 
of the weak and strong solutions to this model. Moreover, we establish the existence of a finite dimen- 
sional global attractor to this dissipativc evolution system. Preliminary computational tests indicate 
that our hyper— diffusion model does not exhibit any of the nonphysical instabilities near the literal 
boundary which are observed numerically in other models. 



AMS Subject Classifications: 35Q35, 65M70, 86A10. 

1. Introduction 

The planetary geostrophic (PG) models, the adiabatic and inviscid form of "thermocline" or "Phillips 
type II" equations, are derived by standard scaling analysis for gyre-scale oceanic motion (see [17], [19], 
[20] and [34]). In their simplest dimensionless f3— plan form, these equations are: 



d xP -fv 2 = (1) 

d y p + M = (2) 

d zP + T = (3) 

V • v + d z w = (4) 

d t T + v ■ VT + wd z T = , (5) 



in the domain ft = {(x,y,z) : (x,y) £ M C M 2 , and z e (— h, 0)}. Here V = (d x ,d y ), v = (uijUa) 
denotes the horizontal velocity field, w the vertical velocity, p is the pressure, T is the temperature, and 
/ = /o + (3y is the Coriolis parameter. A diffusive term, K v d^T, is commonly added to equation (5) as a 
leading order approximation to the effect of microscale turbulent mixing. Thus equation (5) becomes: 

d t T + v ■ VT + wd z T = n v d 2 z T. (6) 

In [28] it is argued, based on physical grounds, that in closed ocean basin and with the no-normal- 
flow boundary conditions, this model can be solved only in restricted domains which are bounded away 
from the lateral boundaries, dM x (— h, 0). Thus, it cannot be utilized in the study of the large-scale 
circulation. Furthermore, it has been pointed out, numerically, by [5] that arbitrarily short linear distur- 
bances (disturbances that are supported at small spatial scales) will grow arbitrarily fast when the flow 
becomes baroclinically unstable. This nonphysical growth at small scales is a signature of mathematical 
ill-posedness of this model near unstable baroclinic mode. Therefore, Samelson and Vallis proposed in 
[28] various simple friction and diffusion schemes to overcome these physical and numerical difficulties. 
In particular, they propose a linear drag (or viscosity) in the horizontal momentum equations and a 
horizontal diffusion in the thermodynamic equation (6). Thus the full dimensionless system becomes: 
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dxP - fv 2 = - 


■evi + v v d\v\ + v h A h v x 


(7) 


9yP + fVl = - 


ev 2 + v v d 2 z v 2 + v h A h v 2 


(8) 


d zP + T = 




(9) 


Vv + d z w = 




(10) 


dtT + v • VT - 


- wd z T = n v d 2 z T + n h A h T , 


(11) 



where A/, = -^i + Jp-- Notice that the incompressibility and hydrostatic balance are retained in the 
above system. Here, the coefficients Vh, v Vl e and Kh are non- negative and small. 

In the case where Vh > 0, v v > 0, e > 0, k v > and Kh > 0, with the ocean being driven by the wind 
stresses at the top surface, and with no-slip boundary conditions and no-heat fluxes on the side walls, 
T s = dM x (— h, 0), and at the bottom, M x {z = —h}, the above system has been studied analytically 
by [25] and [26]. Under this type of boundary conditions the first two authors have improved in [2] 
the results of [25] and [26] and established global regularity and well-posedness to the system (7)-(ll). 
In particular, they show in [2] the existence of finite dimensional global attractor to this viscous three 
dimensional planetary geostrophic (PG) model. This global well-posedness result provides a rigorous 
justification to the scheme suggested by Samelson and Vallis [28] under the above conditions on the 
coefficients. 

On the other hand, the case where Vh — 0, v v — 0, e > 0, k v > and Kh > 0, has been suggested in 
related schemes by other authors [12], [21], [22], [36], and [37]. Equations (7) - (11) become 

d x p - fv 2 = -ev 1 (12) 

d y p + fv 1 = -ev 2 (13) 

d zP + T = (14) 

\7v + d z w = (15) 

d t T + v ■ VT + wd z T = n v d 2 z T + n h A h T . (16) 

In this situation, and by following [28], the frictional-geostrophic relations (12) and (13) can be solved 
locally for the horizontal velocities to give: 

vi = -l(ep x + fp y ), v 2 = i(fp x - ep y ) , (17) 

where 7 = (f 2 + e 2 ) -1 . Assuming we have a nice solution up to the boundary, equations (17) and (14) 
imply: 

d z v 1 = - 1 (eT x + fT y ), 8 z v 2 = 7 (/T x - eT y ) . (18) 

As a result of the above, the no-normal-flow boundary condition on the lateral boundary yields the 
following boundary condition on the temperature: 

e— + f— = on the lateral boundary, T s = dM x (-h, 0) , (19) 
on ' os 

where ^ and ^ denote the normal and right-hand tangential derivative, respectively, so that with 
e = (eni — fn 2 , fm + en 2 ) /(e 2 + f 2 ) 1 ^ 2 , we have the following oblique boundary condition on the tem- 
perature. 

dT 

— = 0, on r s . (20) 

On the other hand, the no heat flux at the lateral boundary gives: 

dT 

—Kh-Q— = on the lateral boundary, dM x (— h, 0) . (21) 
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Based on the above it is argued in [28] that in the presence of rotation, i.e., / ^ 0, and in order to be 
able to satisfy both boundary conditions (19) and (21) one has to add to the thermodynamics equation a 
higher order (biharmonic) horizontal diffusion. Otherwise, the problem (12)-(16) subject to the additional 
boundary conditions (19) and (21) is over determined and hence is ill-posed, which explains the cause 
for the observed numerical instabilities near the lateral boundary. In order to resolve this discrepancy 
Samelson and Vallis propose in [28] the following "horizontal hyper-diffusion" planetary geostrophic 
(PG) model: 

dxP - fv 2 = -ev! (22) 
d y P + M = -ev 2 (23) 
d zP + T = (24) 
V • v + d z w = (25) 
dtT + v • VT + wd z T = k v 8 2 z T + Kh A h T - \A 2 h T , (26) 

subject to the no-normal-flow together with the boundary conditions given in (19) on the lateral boundary. 
However, the no-heat flux boundary condition (21) is replaced by the boundary condition: 

dT dT 

-n h — + XA h — = on the lateral boundary, T s = DM x <-h, 0). (27) 
on on 

It is worth stressing that the differences between the boundary conditions (19) and (21) is due to the 
Coriolis parameter. Therefore, it is natural to require A — > when / — > 0. Since the diffusion operator 
K v d 2 T + n h A h T - \A 2 h T in the model (22)-(26) with the boundary conditions (20) and (27) may not 
positive, the numerical instabilities have been observed near the lateral boundary. 

Motivated by [28] we propose below a variant of the system (34)-(43) which is globally well-posed and 
which possesses a finite dimensional global attractor. 

Let O = M x (-/i,0) G K 3 , where M is a smooth domain in K 2 , or M = (0,1) x (0,1), the PG 
equations with friction and diffusion can be written as: 

Vp + fkxv + V = Q (28) 

d zP + T = (29) 

V • v + d z w = (30) 

d t T + v VT + wd z T = K v T zz -V-q(T), (31) 

where D is the friction or dissipation of momentum and V • q(T) — K V T ZZ is the heat diffusion. 

Naturally, there are two friction schemes: one is the linear drag, i.e., V = ev; the other is the 
conventional eddy viscosity, i.e., V = — e(AhAv + A v d zz v) 1 where Ah and A v are small positive constants. 
The PG model with conventional eddy viscosity and the diffusion q(T) = —K h \7T, has been studied in 
[2], [25], [26]. Here we concentrate on the case of linear drag. As explained above, in a closed basin 
and with a linear drag scheme, it is appropriate that V • q be a fourth-order diffusion. Denoting, from 
now on, by A = A/, = d 2 + d 2 . Our goal in this paper is to study the model with V = ev, and with 
"horizontal" hyper-diffusion 

q(T) = XH\7(V ■ (H T VT)) - K h \7T + »VT ZZ , (32) 

where 

H=() /e (33) 
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and H T is the transpose of H, and A, /j, and e are small positive constants. That is we study the following 
model: 

Vp + fk x v + ev = (34) 

d zP + T = (35) 

V • v + d z w = (36) 

d t T + v • VT + wd z T = Q + K V T ZZ - V • (A#V(V • (H T \7T)) - K h VT + p\7T zz ) , (37) 

where Q is a given heat source. 

Now we impose the appropriate boundary conditions to this friction and hyper-diffusion PG model 
(34)-(37). The natural boundary conditions are no-normal flow condition for the velocity field (v,w) 
and non-flux boundary condition for the temperature T (see, e.g., [18], [27], [28], [29]): 

on T u : w = 0; (38) 
on T b : w = 0; (39) 
on T s : v ■ ft = 0, (40) 

onT u :-K v ^-^a(T-T*); (41) 
dT 

on T b : — - 0; (42) 
oz 

on T s : q(T) ■ ft = 0; (43) 
where r„, r& and T s denote the boundary of f2 defined as: 

r„ = {(x,y,z)er!:z = 0}, (44) 

T b = {(x,y,z)eil:z = -h}, (45) 

T s = {(x,y,z)en:( Xl y)edM}, (46) 

a is a positive constant, ft — (ni,n 2 ) is the normal vector of T s , and T*(x,y) is a typical top surface 
temperature profile. Furthermore, the no-normal flow condition implies a boundary condition on the 
temperature, as explained in (17), (18), and (19): 

dT , , 

^ = 0, onT s , (47) 

where e = (eni — Jn 2 , fn\ + en 2 ) / (e 2 + f 2 ) 1 ^ 2 - Therefore, no-normal flow (47) and no heat-flux bound- 
ary conditions (41)-(43) are natural and proper boundary conditions for the PG model (34)-(37). Since 
(47) is in fact a boundary condition on the temperature, thus, any additional and incompatible boundary 
condition to the temperature, on top of (41)-(43) would make the fourth-order diffusion V • q overdeter- 
mined. Finally, the model is supplemented with the initial condition: 

T{x,y,z,0)=T (x,y,z), (48) 

where T is a given function. 

We observe that in the case when (3 = 0, i.e. / = /o, our "horizontal" hyper-diffusion term V • q(T) 
reduces to the form AA 2 T — p,Ad 2 T— K h AT, which is in the spirit of the hyper-diffusion term in (26) that 
was proposed by Samelson and Vallis in [27] and [28]. However, due to the fact that (3 ^ our proposed 
hyper-diffusion term takes a more involved form, which is necessary to guarantee the dissipativity of this 
operator under the given physical boundary conditions (38)-(43). It is worth adding, that the present 
PG formulation (34)-(43) has been explored with some preliminary computations by Samelson [24] using 
a modified version of the (3— plane numerical code developed in [28]. In the (3— plane case (for which the 
horizontal coordinates are Cartesian x and y, and the Coriolis parameter / = /o + (3y, with /o and f3 
constant), the additional horizontal diffusion terms in the thermodynamic equation reduce to a single 
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term proportional to T xx . The preliminary computations indicate that inclusion of these additional terms 
only slightly modifies the previous numerical solutions reported in [28] which uses hyper-diffusion term of 
the form suggested in (26). With the rigorous analytical results proved here for the modified system(34)- 
(43), this provides new theoretical and mathematical support for the approach and analysis of [27]. Note 
that these preliminary computations did not include the additional mixed horizontal-vertical diffusion 
term /jAT ZZ} which is a crucial term for our rigorous mathematical analysis. 



2. Preliminaries 

It is natural to assume that T* satisfies the compatibility boundary conditions: 

dT* 



= 0, 



on dM. 
on dM. 



(49) 
(50) 



de 

q(T*)-n = 0, 

Let T = T — T*. Due to the compatibility boundary conditions (49) and (50), it is clear that T satisfies 
the following homogeneous boundary conditions: 



dT 

dz 



z=—h 



z=0 



dT 

0; ar 



0; q(T) ■ n 



0. 



By replacing T by T, we have the following new formulation of the system (34)-(48): 
dtf + V • q(f) - K v Tzz +v-VT + wd z f + v-VT* = Q*, 



d ' T+ K? 



0; d z T 



z=0 



= 0; 



dT 

~di 



= 0; q(T)-n 



r. 



= 0, 



(51) 

(52) 
(53) 

(54) 



T(x,y,z,0) =T (x,y,z) -T*(x,y). 

From now on q(T) is given by (32). Here v and w are determined by the use of (18) and (36), and the 
fact that the average of v in the ^-direction is zero. This can be easily seen by integration of (34) and 
(35), which yields the following system: 

Vp + fk x v + ev = 0. (55) 
V-tJ=0, (56) 
v-n = 0. (57) 

Here v and p are the averages of v and p in the ^-direction. Hence, multiplying by v and integrating over 
M, we obtain that v = 0. The exact expressions of v\, vi and w in terms of T are: 



Vl 



i: 



eT x (x,y,Z,t) + fT y (x,y,Z,t) 



e 2 + p 



h\ e T*( Xl y) + fT;(x,y) 

£ 2 + f 2 



V 2 = 



r r) 



-h 



eT x (x, y, t) + fT y (x, y, £, t) 



e 2 + p 



-fT x (x,y,£,t) + eT v {x,y,i,i) 



d£ + [z + 



ff] 



£ 2 + p 



-fT x (x,y,Z,t) + eT y (x,y,Z,t) 



-fT*(x,y)+eT;(x,y) 

e 2 + p 



(58) 



(59) 



e 2 + p 



Using (36), (38) and (39) we have 



w(x,y,z) 



[ {V-v)(x,y, 
J-h 



(60) 
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From (58), (59) and (60) we obtain 



w 



+ 



IT 

j-h j-h 



+ 



+ 



(e 2 + /2) 2 

z + h r° n 

h J-h 



eAT(x,y,£,t) - f T x (x,y,£,t) 

£ 2 + / 2 

(-fT x (x,y,£,t) +eT y (x,y,£,t)^ 



(61) 



d£_drj 



eAT(x, y, £, t) — f T x (x, y, £, t) 



e 2 + f 2 



(e 2 + /2) 2 

z(z + h) 



(-fT x (x,y,£,t) + eT„(s,y,£,t)) 



d^d-q 



eAT*(x,y)-f T*(x,y) 2fof, { 
+ , o. ,n,n {-jT x {x,y) + eT y {x,y)) 



€ 2 + p 



(e 2 + /2)2 



We denote by L P (Q) and H m (Q) the standard LP— spaces and Sobolev spaces, respectively. Following 
the notations in [2] and [25], we set 



V = \ R G C°°(ft) : ( + — T 



a 



= 0; 9 Z T 



z=0 



z——h 



= 0; g(T) • n 



= 



For every R G C°°(fl), denote by 



l^lk - (pl| 2 ffl( o) + l|Vi2||^ 1(n) ) 



and 



\ r \\h^(Q) + \\ AR \\ 2 H2(n) 



1/2 



1/2 



\\R\\v t = ( 

Also, we denote by 

V 2 = the closure of V with respect to the topology induced by the norm || • ||y 2 . 

and 

V4 = the closure of V with respect to the topology induced by the norm || • \\v 4 - 

It is clear that V2 and V4 are separable Hilbert spaces. Next we define the bilinear form a : V2 x V2 
is defined by 



(62) 

(63) 

(64) 

(65) 
(66) 



a(Ri,R 2 ) = a 

Jz=0 



Ri R2 + -p-ViJi • Vi? 2 



dxdy + (67) 

+ / [K h VRi ■ VR 2 + K v (Rx) z (R 2 ) z + AV • (H T VRi) V • {H T VR 2 ) + {iV(Ri) z ■ V(ifc)*] dxdydz. 
We will denote by 



|i?| 2 dxdydz, 



\R\\ 2 = a(R,R) = a 



z=0 



|i?l 2 + £|vi?| 2 



dxdy 



(68) 
(69) 



+ / [K h \VR\ 2 + K V \8 Z R\ 2 + \\V ■ (H T VR)\ 2 + n\VR z \ 2 ] dxdydz. 
Now we give the definition of weak and strong solutions to the model. 
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Definition 1. Let S be any fixed positive time. A function T(x,y, z,t) is called a weak solution of 
(52)-(54) on [0, S] if 



T G C w ( [0, S] , L 2 (fi) ) n L 2 ( [0, S] , V 2 ) , T t G L 1 ( [0, S] , V 2 ') , 



and 



T(t)^ dxdydz- ( T(t )ip dxdydz + [ a(T(s),ip) 
t 



ds + 



(70) 



t Jn 



v-VT(s)+wT z (s)+vVT* 



tp dxdydz ds 



Q*ip dxdydz ds, 



t Jn 



for every ip € Vi. Here C w ([0, S], L 2 (0)) is the functional space of all weakly continuous functions from 
[0, S] to L 2 ( fl). Furthermore, T(x, y, z, t) is a strong solution of (52)-(54) on [0, S] if it is a weak solution 
and 

f(x, y, z, t) G C([0, S],V 2 ) n L 2 ([0, S], V 4 ). 
Next, we give some remarks about the following boundary value problem: 
( V • q(R) - K V R ZZ = g(x, y, z) in ft 



d z R+^R 



o; d z R\ z= _ h 



0; l| rs =0; q(R)-n\ Fs 



0, 



(71) 



where q(R) is given by (32). By integrating by parts and the boundary conditions (51), we have 

J2i(V ■q(R 2 )-K v {R 2 ) zz ) dxdydz = [ R 2 (V • q(R{) + d z q 3 {R{)) dxdydz = a(R u R 2 ), (72) 



for all Ri,R 2 G V. Notice that 

e • n = e ^ 0. 

Namely, the vector e is not tangent to dM. Using the classical regularity results in smooth domains, 
Mcl 2 , for the Laplacian operator with oblique boundary condition, we have 

Proposition 2. Suppose that R(-, ■, z) satisfies the boundary condition ^§| r = and that AR(-, ■, z) is 
in L 2 {M) for every fixed z. Then there exists a constant C > (independent of z) such that 



\\R(;;z)\\ <C\\\AR(;;Z)\\ 2 L2{M) + \\R(;;Z)\\ 

l 2 (m) 

Consequently, there is a constant Co > such that 

^r\\R\\ 2 <\\R\\ 2 v 2 <C \\Rf. 



(73) 



(74) 



Proof. As we mentioned earlier the proof of (73) is a result of classical regularity theory of elliptic 



, OR 
equations. Notice that since — — 

oe 



— (thanks to (71)) we also have 



dR> 



8e 



0. Now applying (73) 



to d z R instead of R with the corresponding boundary condition 



OR, 



de 



= implies (74). 



□ 



Using (74), Lax-Milgram Theorem shows that there is a unique solution R G V 2 for the boundary 
value problem (71) satisfying 

11*11 < \g\- (75) 
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Moreover, using integration by parts in (71) leads to 

| ff | 2 = |V • q(R) - K V R ZZ \ 2 = [ [V • q(R) - K V R ZZ } 2 dxdydz 

Jn 

= I [\V ■ q{R)\ 2 + K 2 V \R ZZ \ 2 + 2K V V ■ q{R)R zz ] dxdydz 
Jn 

= ( [\V ■ q(R)\ 2 + K 2 \R ZZ \ 2 + 2K v K h \VR z \ 2 + 2K V ^\VR ZZ \ 2 + 2K V X\W ■ (H T WR Z )\ 2 ] dxdydz + 
Jn 

+2a [ [ K h \\7R\ 2 + 2A|V • (H T \7R)\ 2 ] dxdy. 
As a result, we have 

\Rzz\<C\g\. (76) 

Notice that the principle part of operator V -q(R) —K V R ZZ is — \A(-AR + R) — K V (—AR + R) zz . Then, 
using a symmetry argument in the z-direction and the standard regularity results for the Laplacian 
operator (see, for example, [13] p. 89), we get 



AR + R\\ HHn) <C\g\. 



Therefore 



\R\ 



(l|Ai?|| 2 ff2(0) 



\ 1/2 

) <C\g\. 



(77) 



As a result of Proposition 2, (74) and Rellich Lemma [1], one can show that the operator V-q(R) — K v R zz 
with domain V4 is a positive self-adjoint operator with compact inverse. Therefore, the space L 2 (Vl) 
possesses an orthonormal basis {4>k{x, y, z)}j° =1 of eigenfunctions of the operator V • q(R) — K V R ZZ , such 
that 



V • q{(j>k) - K v ((f> k ) zz = X k 4>k, 



(78) 



where < Ai < A 2 < • • • , and lim A& = 00. Moreover, by stand results(cf. for example, [4]), we have 

k^oo 



Ci " Ai ' 



(79) 



We will denote by H m = span{0i, • 
onto H m . 



and by P m : L 2 (fl) — ► H m the L 2 (Q) orthogonal projection 



For convenience we recall the following classical inequality about the trace operator: 



\R\ 



< C2\\R\\ H i(a), 



(80) 



and the following version of Sobolev embedding and interpolation Theorems (cf. for example, [1]): 



\\h(x,y)\\ LHM) < C3||/.(x,y)||^ 2 (M) ||^,2/)||^(M)' yhE H ^ M ^ 
\\K x ^v)\\lo(m) < C q \\h(x,y)\\ H i( M) , V h £ H^{M), and q < oo, 

\\h(x,y)\\ L ~ (M) < C3||/.(x,2;)||^ 2 (M) ||Mx,y)||^ 2 (M) , Vfte H 2 (M), 

\\h{x,y)\\ m(M) <C3||M^,y)l!L 2 (M)ll^(^y)ll^ 2 (M)' VheH 2 (M), 

\\Ah(x,y)\\ LHM) ^CsllM^yJII^MjIIM^tf)!!^)' VheH 2 {M). 

\\g(x,y,z)\\ L 3 {n) < d\g(x,y, z)\^ 2 \\g(x iy , z)^^, y ^ g H i( Q y 
\\g(x,y,z)\\ L 6 (n) < C 4 \\g(x,y,z)\\ H i {n) , 



1/2 



(81) 



(82) 
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Also, the following integral version of Minkowsky inequality for the L p spaces p > 1. Let fii C K™ 1 and 
f^2 C i m! be two measurable sets, where mi and m,2 are two positive integers. Suppose that f(£,r]) is 
measurable over fii x fi 2 - Then, 



/ (7 /P </ (7 i/(^)r^) /P dv. 

J si! \Jn 2 / J -^n 2 \Js2i / 



(83) 



Hereafter, C, which may depend on the domain O and the constant parameters e, /o, /3, a, if/,, lf„, A, /i in 
the system (52)-(54), will denote a constant that may change from line to line. 



3. Global Existence, Uniqueness and well-posedness of Weak Solutions 

Now we are ready to show the global existence and uniqueness of weak solutions to the system (52)- 
(54). 

Theorem 3. Suppose that T* GjP(M) and Q e L 2 (fl). Then for every T (x, y, z) G L 2 (ti), and S > 0, 
there is a unique weak solution T of the system (52)-(54) . Moreover, T satisfies 



d t TeL^(0,S;Vl), Vr>l, \T(a)\ 2 + f \\T(s)\\ 2 ds < K 2 (S,Q,T ,T*), Vae[0,5], (84) 

Jo 



where V{ is the dual space of Vi, K\(S,Q,To,T*), and Kz(S,Q,To,T*) are as specified in (100) and 
(102), respectively. 

Proof. First, let us prove the existence of the weak solution for system (52)-(54). We will use a Galcrkin 
like procedure, based on the eigenfunctions {4>k} kLi, to show the existence. Let m £ Z + be fixed, the 
Galcrkin approximating system of order m that we use for (52)-(54) reads: 



d_ 

dt 



+ V-q(T m )-K v (T m ) zz + P n 



v-VT m + w^ + {v-V)T* 
oz 



T m {x,y,z,0) = P m [T Q (x,y,z)-T*(x,y)}, 



(85) 
(86) 
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where T m = X^fcLi a k(t)(f>k(x, y, z), and v(T m ) = (i>i, v 2 ), w = w(T m ) given in terms of T m by the formulas 
below, and 



q(T m ) = -K h VT m + AHV(V • (H T VT m )) + M V(T„ 



e{T m ) x {x, y, i, t) + f(T m ) y {x, y, £, t) 



e 2 + P 



Z+ 2 



h\ eT*(x,y) + fT*(x,y) 



e 2 + f 2 



r r) 



e(T m ) x (x, y, £, t) + f{T m ) y {x, y, £, t) 



'"2 



-f 

J-h 



h J-h ^+P 
f(T m )x(x,y,£,t) + e(T m ) y (x,y,t,t) 



d£+[z + 



r r) 



e 2 + f 2 

-f(T m ) x (x, y, £, t) + e{f m ) y (x, y, £, t) 



h\ -fT*(x,y) + eT;(x,y) 
e 2 + / 2 



(87) 
(88) 

(89) 



h J-h 



e 2 + p 



d£dr/, 



w = 



z r r) 



+ 



hJ-h 
2/o/ 



eAT m {x, y, £, t) - f Q {T m ) x {x, y, £, t) 



+ 



e 2 + p 

, , r2 , 2 (-f{T m ) x (x,y,U) + <T m ) v (x,y,U)) 



(90) 



(e 2 + P) 

z + h r° n 



d£dr] + 



2/o.f 



eAT m (x, y, £, t) - f (T m ) x (x, y, £, t) 



(£ 2 + /2) 2 

z(z + h) 



c 2 + P 

f(T m ) x {x,y,€,t) + e(T m ) y (x,y,€,t) 
eAT*(x, y) — f T*(x, y) , 2/ / 



d^drj 



e 2 + P 



+ 



(e 2 + /2)2 



{-fT*(x,y) + eT;(x,y)) 



(91) 



We stress again that v and w depend on m since they are functions of T m . However, we will drop the 
explicit dependence m to simplify the notation. The equation (85) is an ODE system with the unknown 
afc(t), k = 1, ■ • • , m. Furthermore, it is easy to check that each term of equation (85) is locally Lipschitz 
in T m . Therefore, there is a unique solution ak{t), k = 1, ■ • • , m, to the equation (85) for a short interval 
of time [0, S*). By taking the L 2 (f2) inner product of equation (85) with T m , we obtain 



ld\T„ 



2 dt 



p 



n 



v ■ VT m wd z T m 



(v ■ VT*)T m dxdydz + 



(92) 



T m dxdydz 



[ T m Q* 



dxdydz. 



It is easy to show by integrating by parts and by using the boundary conditions (38)-(40) that 



/ 

Jn 



v • VT m + wd z T„ 



T m dxdydz = 0. 



(93) 



Furthermore, by Holder inequality we have 

/ (v ■ VT*)T m dxdydz 



<C||w|U6 ( n)||VT*|| L 3 ( n) \T m \. 



By (82), (88) and (89), we have 

IMU°(fi) < C\\v\\ H i (n) < C 



+ \\T*\ 



H 2 (M) 



< c 



+ \\T*\ 



H 2 (M) 



(94) 
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By (81), we obtain 

||VT*|| L 3( M ) < C||T*|| ff 2( M ). 
As a result ol the above estimate, we get 

(v ■ VT*)f m dxdydz 



< C 

<I 
~ 4 



T 



IT* I 



H 2 (M) 



+ 



||T*||_H- 2 (M) |T TO |. 
C||T*|| 2 ff2(M) (l + |T m | 2 ). 



(95) 



Applying Cauchy-Schwarz inequality and the definition of Q*, we obtain 

Q*T m dxdydz 



< 



J 

Jn 



QT m dxdydz 



+ 



[ 

Jn 



V • q(T*) T m dxdydz 



< \Q\ |T m | + C||T*|| H 2( M )||T TO || 

<\\Q\ 2 + \\f m \ 2 + c\\T*r HHM) + \\\f„ 

Therefore, from the above estimates (93)-(96) and (74), (92) gives 
d\T m f 



dt 



< C 



|Q| 2 + ||T*" 2 



H 2 (M) 



C\\T*\\ 



H 2 (M)\ ± m\ 



Thanks to Gronwall inequality, we conclude 



\T m (t)\ 2 < [|T | 2 + C(|Q| 2 



+ \\T* 



C\\T 



* II 2 



\H 2 (M) 



(96) 



(97) 



(98) 



\H 2 (M) j 

when < t < S*. But since the right hand side is bounded as t goes to S*, we conclude that T m {t) must 
exist globally, i.e., S* = +oo. Therefore, for any given S > and any t e [0, S], we have 

|T m (t)| 2 KK^Q^T*), 

where 



|T | 2 + C 



(\Q\ 2 



\T*\\h 2 (m) 



tfi(S,Q,T ,T*) 

By integrating (97) with respect to t over [0, S], and by (99), we get 

|T m || 2 d S <if 2 (5,Q,T ,T*), 



C||T*|| 

H 2 (M) 



(99) 
(100) 



./o 



where 



K 2 (S,Q,T ,T*) = \T \ 2 + C IQl' + llT^^ + llT^^K^Q^TolT*) 



(101) 



(102) 



and K\(S, Q,Tq,T*) is as in (100). Notice that the estimate (99) is unbounded in time (i.e., as S — ► oo), 
but it is uniformly bounded in to. However, in Section 5 we will present a sharper estimate which is 
asymptotically bounded in time. As a result of all the above we have T m exists globally in time and is 
uniformly bounded, in to, in L°°([0, S]; L 2 (f2)) and L 2 ([0, S]; V2) norms. 

Next, let us show that dtT m is uniformly bounded, in to, in the L 3 ^- 1 ([0, S 1 ]; V 2 ') norm for every r > 1. 
From (85), we have, for every ip € V2 



~ r)T 
v-VT m + w°-^- + (v ■ V)T* 
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Here, (•, •} is the dual action of V{. It is clear that 



\(P m Q*^)\ < \{P m QM + |(V • q(T*)^)\ 

< \\QM + C\\T*\\ HHM) Ml, 
and by integration by parts we have 

[-K v (T m ) xx + V-q(T m ),1>)\<C\\T, 
Next, let us get an estimate for 

df, 



v-V(T m + T*)+w 



dz 



v-V{T m + T*)+w 



dz 



ip dxdydz 



v-V(T m + T*) + w 



dT m 
dz 



ip m dxdydz 



where ip m = P m tp- Thus, by integration by parts, we obtain 

df 



P 

± rr 



v ■ V(T m + T*) + w- 



dz 



V ■ \7lpm + W 



dlpm 

dz 



(Tm + T*) dxdydz 



Next, we estimate 



Jjv-V^m)(Tm + T*^ dxdydz 



Applying (63), (82) and Proposition 2, we have 

||VV>m||l,6(fi) < CllVV'mlljji(n) < C||^ TO ||. 

By (82) and (94), we reach 

||«IUs ( n) < C\\v\\^ m M 1 / 2 < C \\\fm\\ + ||T*|| ff2(M) 



Since 



+ T*\ < \T m \ + h 1 ' 2 \\T*\\ 



the above inequalities imply 

f (V ■ V^m) (' 



< c 



T m + T* 

\\T m \\ + \\T*\\h 2 (M) \T m \ + \\T*\\ L 2r M ) \\i>m\\- 
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For the other term in (106), we use (60) to get 
J wd z ^> m (r m + T*) dxdydz 

J^j ' V • t> de) d z yj m (f m + T*) dxdydz 
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< 



M 



|V • v(x,y,£,t)\ d£ ) \d x ip m (x,y,z)\ T m (x,y,z,t) +T*(x,y) 



dxdydz 



-h 



\V -v(x,y,t,t)\ d£ / \d z ip m (x,y,z)\ T m (x,y,z,t) +T*(x,y) 



-h 



dz 



dxdy. 



Using Cauchy-Schwarz inequality, we obtain 



/ \d z ip m (x,y,z)\ T m (x,y,z,t) + T*(x,y) 
J-h 



dz 



< 



/ \d z ^ m {x,y, 
J-h 



z)\ dz 



T m (x,y,z,t) + T*(x,y) 



dz 



Thus, 



wd z tp m (r m + T*^j dxdydz 



< 



< 



M 



\V-v{x,y,Z,t)\d£ 



\d z ip m {x,y,z)\ 2 dz 



-h 



\T m {x,y,z,t) + T*(x,y)\ 2 dz) dxdy 



'us: 



|V • v(x,y,£,t)\ dn dxdy 



M \J-h 



\d z ip m (x,y,z)\ dz) dxdy) x 



M 



T m (x,y,z,t) + T*(x,y) 



dz ) dxdy 



(108) 



here we apply Holder inequality and 1/r + 1/r' = 1, and r > 1. By using Minkowsky inequality (83), we 
get 



M 



(x,y,z,t) + T*{x,y)\ 2 dz) dxdy 



< 



f ( f \T m (x,y,z,t)+T*(x,y)\ 2r dxdy) dz 

yJ-h \JM / j 

Thanks to (81), for every fixed z and t we have 



Jm 



x,y,z,t)+T*(x,y) 



2r' 



)2r 
< c r 



T m (z,t)+T* 



L 2 (M) 



T m (z,t)+T* 



H 2 (M) 
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As a result of the above and Holder inequality, we obtain 



o 

-h \JM 
■0 



T m {x,y,z,t) + T* 



2r 



dxdy I dz 



< C r 



< C r 



-h 



T m (z,t)+T* 



L 2 (M) 



dz 



T m (z,t) + T* 



H 2 (M) 



dz 



T m + T* 



T m + T 



v 2 



(109) 



Similarly, by using Minkowsky inequality (83) and (81), we get 

(/ (/ \ dz ^ m ( x,y,z ^ ' dz j dxd y^ 

~ (/ (/ \ dz ^ m ^ x ' y,z ^ 2r dxd y^j d ^j 

( f° 2 A' 7 ' 

~ C \J h " a ^" l "« 1 ( M ) 
<C r U m \\. 

Here we used Proposition 2. As a result of (108), (88), (89), (109) and the above estimate, we have 

/ wd x il> m (f m + T*) dxdydz (110) 
Jo. 

r+l 3r— 1 

< C r (|T TO | + ||T*|| £ 2( M )^ (l|r m || + \\T* \\ H 2 (M)^j \\4'm\\- 

By (107) and (110), we have 

df. 



Prr 



v-V(T m + T*)+w 



dz 

r+l 



<a(|r m | + ||T*|| i2(M) ) 2r (|jr m || + ||r*|| ff2(M) ) 2r 

Since ip € V2, then the Fourier series 



V] ( / V^fc dxdydz \ (f) k = ip m + V" ( / V^fc dxdydz ) <£ fc 

fc=l ' fcm+l / 



converges to ^ m V2 (cf. [13] p. 64). As a result, we get, 

||VUv 2 <c\\n V2 . 

Therefore, from the above and Proposition 2 we have 



v ■ V(T TO + T*) + w- 



dz 



(111) 



r+l . 

< C r (\T m \ + ||T*||i2( M )^ (||r m || + ||T*||#2(M)) 
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for every r > 1. By the estimates (104)-(111), (99) and (101), we have 



d t T m ,i> 



<C(\Q\ + \\T*\\ mM) )U\\+C\\T, 



+a( 



T m \ + \\T* 



\L 2 (M) 



+ \\T* 



Thus, due to (99) and (101), we have 



f 

Jo 



IdtT^t)^,- 1 dt < K 3 (S,Q,T ,T*,r), 



where 



K 3 (S,Q,T Q ,T*,r) =C[l + \Q\ + \\T*\\ H * {M) ] 3 - 1 S + CK 2 (S,Q, \T \,T*) + 
+C, 



(112) 



(113) 



Ki(S,Q,T ,T*) + ||T*|| ff 2(Q) 



K 2 (S,Q,T ,T*) + \\T*\\jp 



*„2 g 



(M) 



Therefore, dfT m is uniformly bounded, in to, in the L 37 — 1 ([0, S]; V 2 ) norm, for every r > 1. Thanks to 
(99), (101) and (112), one can apply the Aubin's compactness Theorem (cf., for example, [7], [15], [31]) 
and extract a subsequence {T mj } of {T m }; a subsequence {v mj } of {v m — v(T m )} and a subsequence 

{d t f mj } of {dtf m }; which converge to f e i°°([0, S]; L 2 (n)) n L 2 ([0, S]; V 2 ) and <9 t T 6 ([0, 5]; V£), 

respectively, in the following sense: 

" T m . -T inL°°([0,S];L 2 (fi)); 
T m] ^T in L 2 ([0,S};H 1 (n)) weakly; 
a t f m . -» ftT inl^ ([0, 5] ; KO weakly. 

Notice that since T m . e V , by integration by parts it is clear that 

/ f mj (x,y, z,t)P mj ip dxdydz - / f m] (x, y, z, t )P mj tp dxdydz + a(f mj ., P mj tp) 
Jn Jq 



+ J t J n [( v ™ 3 ■ VT m .) P m ^ + w mj d z T mj P mj ip + (v mj ■ VT*)P„ 



,1> 



> t 

= / Q*P mj ip dxdydz, 



dxdydz 



for every ip <E C°°([0, S]; V 2 ), and for every t, to s [0, £]. By passing to the limit, one can show as in the 
case of Navier-Stokes equations (see, for example, [7], [31]) that T also satisfies (70). In other words, T 
is a weak solution of the system (52)-(54). 

Next, we show the uniqueness of the week solution. Let T\ and T 2 be two weak solutions of the system 
(52)-(54) with respect to initial values % and f$, and let (v' = (vi,v 2 ),u/) and (v" = (v'{,v%),w") be 
given by (58)-(61) with respect to Ti and T 2 . Denote by u = v' — v", u 3 = w' — w" and X = T\ — T 2 . 
Unfortunately, since we only have 

d tX eL^(0,S;Vl), Vr>l, 

we are not able to apply the standard energy method and Lions Lemma (cf. [31] Lemma 1.2. p. 260). 
However, let us instead consider 



Xm = PmX, and Xm = X - PmX 
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where P m is the L 2 (£l) orthogonal projection onto H m . It is clear that %m satisfies the following equations: 



dtXm + V • q( Xm ) - K v { Xm )zz + P m [u ■ VTi + u 3 0*Ti + v" ■ Vx + w" ' Xx + u ■ VT* 
Xm(x, y, z, 0) = P m (f^(x, y, z) - f^(x, y, z)j , 

where u = (ui,u 2 ) and 



q(Xm) = -K h V Xm + XHX7(V . (H T \7 Xm )) + l^{Xm)z 
z £Xx(x, y, g, t) + f Xy (x, y, g, t) 
h 



Ul 



(x,y,z,t) = I 



e 2 + p 

h J-hJ-h 



d^dz, 



u 2 (x,y,z,t) 

r0 



£ 2 +/ 2 

/Xz(z, V, £, t) + e Xy (x, y, £, i) 



+ / 2 

1 A" A 2 -fXxjx , y, g, t) + ex^, y, g, t) 



d£dz, 



u 3 (x,y,z,t) = 
Vof 



Z rT) 

-h J-h 



e 2 +f 2 

eA X (x, y, £, i) - /oXx(z, y, £, *) 



+ 



£ 2 + /2 

{-fXx{x, y, 1 1) + e Xv {x, y, £, i)) 



+ 



(e 2 + / 2 ) 2 

z + h f° f z UAx{x,y,€,i)- foXx(x,y,£,t) 



d^d-q + 



-hJ-h 



e 2 + / 2 



2/o/ 



:{-fXx{x,y,£,t) + exy{x,y,£,t)) 



(e 2 + /2) 2 

By taking the V^' dual action to equation (114) with Xm> we obtain 
MxJ 2 



d£dz. 



2 efe 



+ llXr. 



u • VTi + u 3 a,Ti + v" • V X + w"x, + u • VT*, Xr, 



u • VTi + u 3 fl z Ti + w" • Vx + u/'x, +«'Vr 



Xm dxdydz. 



||Vr*|| L 6 ( n)||«|U3( n) |Xr. 



Next we estimate the equation (120) term by term, 
(i) 

/ u ■ VT* Xm dxdydz < 
Jn 

From (117) and (118), by applying (82) and Proposition 2, we have 

ll«IUs (n) < C||Vx|ks (n) <C|| X ||. 
Thus, by the above and Proposition 2, we obtain 



u ■ \7T* Xm dxdydz 



<C||T*|| ff2(M) (| Xm | HxJI + IXml \\Xr, 



(114) 
(115) 

(116) 
(117) 



(118) 



(119) 



(120) 



(121) 
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L 



u ■ VTiXm dxdydz 



< C 



In (/ l Vx ( x ' y ' z '^l dz ) \VTi\fay>z>t)\x m (x,y,z,t)\ dxdydz 



M 



Wx{x,y,z,t)\dz ) / |VTi(x,y, |xm(x,y, z,t)\dz ) dxdy 



Applying Cauchy-Schwartz inequality to the above, we obtain 



/ 

Jn 



u ■ WTiXm dxdydz 



~ C Im (/j Vx|dZ ) (/j VTl|2dZ ) 2 (/j Xm|2dz 



dxdy. 



By using Holder inequality we reach 



/ 

Jn 



u ■ VTiXm dxdydz 



< C 



M 



I if 

JM \J-h 



iVTxfdz) dxdy 



f IVxh 
J-h 

I (f 

JM \J-h 



dz I dxdy 



dz I dxdy 



By using Minkowsky inequality (83) we get 

I i^J \Xm{x,y,z,t)\ 2 dz^j dxdy 
Thanks to (81) for every fixed z we have 



< 



h \JM 



\x m {x,y,z,t)\ dxdy dz. 



M 



\Xm(x,y,z,t)\ dxdy] < C \\Xm(z, t) ||£ 2(M) \\Xm(z, t)\\ 2 H 



<C\\ Xm {z,t)\\h { M) \\Xm{z,t)\\m { M)- 

As for estimate (73), it is easy to obtain 

< C (\\A Xm (z,t)\\L'(M) + \\Xm(z,t)\\L*(M)) ■ 

As a result of the above we reach 



4 \ 4 a 1 

\ X m(x,y,z,i)\ dxdy ) < C \\ Xm {z,i)\\l 2{M) \\Ax m (z,t)\\* 2(M) + C\\x m {z, 



M 



lL 2 (M) 



li 2 (M)' 



Applying Young's inequality, we obtain 



h \JM 



\Xm{x,y,z,t)\ A dxdy ] dz 



<C I \\Xm(z,t)\\l 2{M) \\AXm(z,t)\\l 2(M) dz+\ Xm (z,t)\ 



L 2 (M) ' 




< C 



h \\Xm{z,t)\\ 2 L2{M) dz^j Q \\Ax m (z,t)f L2{M) dz 



<c\ Xm (tT*\A Xm (t)\-* +c\ Xm (t)\ 2 - 
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Therefore, 



/ (/ \ Xm ( x,y,z,t ^' 2 dz ) dxdy 



< C \ Xm (t)\ 3/4 |A Xm (i)| 1/4 + C\ Xm (t)\ 2 . (122) 



As for obtaining (122), by using (83), (81) and (73), and proposition 2, we have 



/ (7 \V X (x,y,tt)\dC\ d.nlij 

JM \J-h 



As a result of (122) and (123) we have 



<C|x| 1/4 |Axr+C|x| ,/2 |lxll' /2 . (123) 



/ 

J 9. 



u ■ VTiXm dxdydz 



< C 



Ti 



(iii) Notice that from (60) we have 

By integrating by parts, we get 
"3 d z fi Xm dxdydz 



(\Xm\ \\Xm\\ + \Xm\ lfi \\Xmf /4 \ X m 1 3/4 1 1 Xm W'") ■ (124) 



= -V- / u(x,y,£,t) d£. 
J-h 



u(x, y, £, t) d^j ■ (xm Vd z T x + dzT^Xm) dxdydz 



By integration by parts we reach 

/ / u(x,y,£,t) d£ ■ Vxmdzfx dxdydz 
Jn J-h 



u(x,y,z,t) ■ Vxr, 



Therefore, 



U3 d z Ti Xm dxdydz 

[ \u\ \V Xm \ \Ti\ dxdyd 
Jn 



u(x, y, £, t)d4) ■ V9 z Xm Ti dxdydz 



+ 



dz I |V9 z Ti \x m \ dxdydz 



u\ dz |Ti| |Va z Xm| dxdydz 



Following the steps of getting estimate (124) we have 



u\ dz ) |V5 z Ti| \x m \ dxdydz 



< 



< C 



J M [J h \ u \ dz ) (/ |V0*Ti| |xm|<fe) dxdy 
(ivf(/j U ' rfZ ) dXdV ) |Vaz?1 ' (j M (/j Xm ' 2dZ ) dxdV ) 

(\Xm\ \\Xm\\ + \Xm\ 1/4 \\Xm\\ 3/i \ X m 1 3/4 || Xm || 1/4 ) , 



\ 1/4 

) 



(125) 
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and 



|U| 

n 

< C 



< C 



|VXm| | Ti| dxdydz 



|V X | dz 



M \J-h 




I Ti| |Vxm| dz dxdy 



M 



|V X | dz max |Ti(-,-,2,t)| / |VXm|^ 

-/i<z<0 



< 



C / T |Vx| dz f /° |Ti| 2 dz) ( [ |(Ti) z | 2 dz) f |V Xm | dz 

JM J-h \J-h / / J-h 

lx™l 1/2 llx m || 3/2 + lx m | 1/4 ||x m || 3/4 lxJ 1/4 llx m || 3/4 



dxdy 

1/4 -o 



dxdy 



Ti 



Ti 



< C 
and 

/ / M dz |Va,Xm| |Ti| dzdydz 

< c [ f |v x | dz f /° |va zXm | 2 <fe) f f |Ti| 2 ^ 

JM \J-h J \J-h 

f (f \V X \dz) dxdy) ( I ( f \T t \ 
Jm \J-h J / \ijif \J-h 



1/2' 



<c\vd zXr , 



s 2 

2 



dxdy 
dz I dxdy 



1/4 



< c 



Ti 



Ti 



4 IxNIxll 1 llXr, 



(iv) By integrating by parts and (54), we reach 

[v" ■ Vx + w"d z x] Xm dxdydz = \ [v" ■ Vxm + w"d z x m ] Xm dxdydz. 

Jq 

Similar to get estimates (124)-(127), we get 

/ [v" ■ Vx™ + w"d z x m ] x m dxdydz 
Jn 

< c [\\xm\\ \f 2 nf 2 \\^ ix m i 3/4 iix m ir /4 + iixmii m^w?^ ix^nixji 174 

+ lx m |^llx m ||^|T 2 r/ 4 ||T 2 || 3 / 4 |x m r /4 ||x m || 3/4 + lx m |'||x m ||'|T2|^||T 2 ||f || Xro ||' 
Therefore, by (121)— (128), and Young's inequality, we have 

d\Xn 12 



dt 
+C 



■<c\\xm\r + 

1 + \\T*\\ 2 HHM) + l|Ti|| 2 + |Ti| 6 ||fi|| 2 + |f 2 | 2 / 3 ||f 2 || 2 l |xJ 2 + 



+C \l + ||Ti|| 2 + |Ti| 6 ||Ti|| 2 + |T 2 | 2 / 3 ||T 2 || 2 1 |x 
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Thanks to Gronwall inequality, we get 



\Xm(t)\ 2 < 



1+ max |T!( S )| 6 ) /Vi( 



max \T 2 (s)\ 



t <s<t 



C 



xe 



xe 



Xm (to)\ 2 + C Axmtoll 2 ds + C 

Jt 

2 / 3 f\\f 2 {s)\\ 2 ds 

Jtg 

max (jT^ + |f 2 ( S )| 2 /3) (/ t * HfxWII 2 + / t * ||f 2 ( S )|| 2 d«) 
) (* - ^) 



'•)|| 2 ds + 



max |xm(s)r 

t a <s<t 



C(l+\\T*\\ 2 H2{Q) 



Notice that the right hand side is bounded, uniformly in m, for every fixed t. By passing the limit, using 
the Lebesgue dominant convergence Theorem, we obtain 



C 



1 + max 

t <s<t 



\x(t)\ 2 < k(to)\ 2 e 

c(l+\\T*\\% Hn) )(t-t ) 



{[T^f + \T 2 (s)\ 2 ^) ||7i( S )|| 2 ds + Jl ||T 2 ( S )|| 2 ds 



xe 



(129) 



Recall that T 1: T 2 e L°°([0, S]; L 2 (n)) nL 2 ([0, 5]; V2), therefore, (129) implies the continuous dependence 
of weak solutions on initial data as well as their uniqueness. □ 

Corollary 4. The weak solution of the system (52)-(54) depends continuously on the initial data. That 
is, the problem is well-posed. 



4. Global Existence, Uniqueness and Well-posedness of Strong Solutions 

In previous section we have proved the existence, uniqueness and well-posedness of the weak solution 
for the reformulated system (52)-(54). In this section we show the global existence, uniqueness and 
well-posedness of strong solutions for the system (52)-(54). 

Theorem 5. Suppose that T* G H A (M) and Q G L 2 (tt). Then for every T (x,y,z) £ V 2 and S > 0, 
there is a unique strong solution T of the system (52)-(54). Moreover, T satisfies 



\\T(t)\\ 2 + f V ■q(T)-K v {T) zz 2 ds<K s (S,Q,T ,T 
Jo 



(130) 



where < t < S and K S (S, Q, Tq, T*) will be specified as in (146). 



Remark: The steps of the following proof are formal in the sense that they can be made more rigorous 
by proving the corresponding estimates first for the Galerkin approximation system based on the eigen- 
functions of operator V -q — K V (T) ZZ with the boundary conditions (53). Then the estimates for the exact 
solution can be established by passing to the limit in the Galerkin procedure by using the appropriate 
"Compactness Theorems" . 

Proof. Let T be the weak solution with the initial datum Tq. we will show that T is a strong solution if 



To(x,y, z) G V 2 . First, we get a priori estimate for 



. Notice that 



- / f zz T dxdydz = / T 2 dxdydz + / T 2 dxdy, 
In Jn Kv J z =o 
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and 



- T, 



11 

aX 



V • q(T) - K V {T) ZZ dxdydz 
K h \VT z \ 2 + K v \f zz \ 2 + A|V ■ {H T VT Z )\ 2 + n\VT zz \ 2 ] dxdydz + 

|V • (T/ T VT)| 2 dxdy + ^ / \VT\ 2 dxdy. 

z=0 Kv 



By taking the inner product of equation (52) with T zz in L (Q), we reach 



z=0 
2/ 



1 d / ~ 2 a 



T 2 cfcciiy 



z=0 



a 



A|V-(7J T VT)| 2 +^|VT| 



+K h \VT z \ 2 + K V \T ZZ \ 2 + A|V • (ff T VT z )| 2 + M |VT Z 



Q* - u • VT - w<9 z T - u • VT* 



T zz dxdydz. 



Let us consider the above equation term by term 
(i) It is clear that 

Q*T ZZ dxdydz 



(Q-V -q(T*))T zz dxdydz 



<C(\Q\ + \\T*\\ HHM) )\T ZZ \. 
(ii) Integrating by parts we get 



-v ■ VT - wT, - v • VT* 



T zz dxdydz 



v z ■ VT + w z <9 z T + w z • VT* + v ■ VT Z + wT zz 

+ 17- I \v -v(t + T*] \ T dxdy. 
J<v J z = L V / J 

Moreover, by (36) and boundary conditions (38)-(40), we reach 

J ■ VT Z + wT zz ^j T z dxdydz = 0. 
By using (58)-(61) and boundary conditions (53), we obtain 



T z dxdydz 



Jn 



v z ■ VT + w z T z + v z • VT* 



T z dxdydz + [ \v ■ V (T + T*)] T dxdy 
K v J z=0 L V /J 



£ 2 + p 



a 

A t> Jz=0 



VT + VT*| 2 T z dxdydz - / V • v \d z T\ 2 dxdydz + 



-h + P 
Thanks to (82) and (73) we have 

e 



1 + J) \X7T(x,y,Z,t)+VT*(x iy )\ 2 dC 



T dxdy. 



VT + VT*| 2 T z dxdydz 



( 2 + f 2 



<C(||VT|| 2 4(0) + ||VT*|| 2 4(M) ^ 



<c(|AT| 2 + ||T*|| 2 ff4(M) ) 
<c(||ff+||T*||^ w ) \f z \. 
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Similarly, thanks to (81) and (73) we have 



1 + J) \VT(x,y,t 1 t) + \7T*(x 1 y)\ 2 d£ 



T dxdy 



K v J z=0 [J_ h e 2 + P 

< C (||Vf||| 4(a) + 

ll^*lli 4 (M)) — o)IU 2 (m) 

< C (|Af | 2 + ||T*||| 4(M) ) ||f (z = 0)|| L2(M) 
<c(||T|| 2 + ||T*|| 2 ff4(M) ) ||r(z = 0)|| L2(M ). 

Following the steps to get the estimate (126) we obtain 

V • v \d z f\ 2 dxdydz 

< C (\AT\ + ||AT*|| L2(M) ) \\f z \\ L2(M) \\VT Z \\ L2{M) 

< C (||f || + ||T*||„ 4(M) ) \\T z \\l* { m) ||Vf z || L2(M) . 
Therefore, from (131)— (136), we have 

|T Z | 2 + J^T 2 dxdy^j + K h \VT z \ 2 + K v \f zz \ 2 + A|V • (H T VT Z )\ 2 + M |VT^| 2 + 
\\V ■ (H T VT)\ 2 + K h \\7T\ 2 ] dxdy 

< C (\Q\ + ||T*||„4 (M) ) |T„| + C (||T || 2 + ||T*|| 2 H4(M) ) |T Z | + 
+C (||T|| 2 + ||T*||| 4(M) ) \\f(z = 0)\\ LHM) + C (||f || + \\T*\\ H * [M) ) ||f,|| L2(M) ||VT,|| L2(M) . 
By Cauchy-Schwartz inequality, we obtain 



(135) 



(136) 



ld_ 

Ydt 



+ 



K v 



d 
It 



T, 



K v 



z=0 



T 2 dxdy) + K h \VT z \ 2 + K V \T Z 



A|V • {H T VT Z )\ 2 + M |VT 2 



(137) 



+TT f U|V-(^ T VT)| 2 + ^|VT| 2 

J^v J z =0 L 



dxdy 



< C (|0| 2 + ||T*|| 2 ff4(M) + ||T|| 2 ) + C (||T|| 2 + \\T*\\ 2 HHM) ) (\\T(z = 0)|| 2 2(M) + |T,| 2 ) . 
Again, by Gronwall inequality, we get 



\T z (t)\ 2 + ^- f T 2 (t)dxdy<K z (t,Q,To,T*), 



where 



K z {t,Q,T ,T*) = e 

x 



C(l + ||T 



*||2 



\h 2 (m)J t + Ki(t,Q,To,T*) 

'roii 2 + iir*ii 2 ff2(M) + ^(t.Q.ro.n + (iqi 2 + \\t*\\ 2 hhm) ) t 

and Ki(t,Q,T ,T*) is as in (100). Finally, let us show that 

f (x, y, z, t) e L°°([0, S],V 2 ) n L 2 ([0, S],V 4 ), 



(138) 



(139) 
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By taking the inner product of equation (52) with V • q(T) — K V T ZZ in L 2 (Q) to get 
M|T|| 2 



2 dt 



+ 



V • q(T) - K V T 2 
w-Vf- wd z f - v ■ VT* 



V • q(T) - K V T ZZ dxdydz. 



Note that 
(i) 



V • q(T) - K V T Z 



dxdydz 



< |Q*| \V ■ q(T) - K V T ZZ \ 

< [IQI + C\\T*\\ HHM) ] |V • q(f) - K v f zz \. 



(ii) 



v ■ VT* 



V • q(T) - K V T Z , 



dxdydz 



<C|«|||VT*|| L co (M) V • q(T) - K V T Z 

<C\\f\\ \\T*\\ HHM) \v-q(f)-K v f Z ; 
(iii) Following similar steps to those which led to the estimate (126) we have 



/ 

Jn 



v ■ VT 



V • q(T) - K V T ZZ dxdydz 



< 



C (|VT| + ||VT*|| L2(M) ) (||AT|| L 2 (a) + ||T*|| ff2(M) ) |v • q(T) - K V T Z 



<C \\T 



ITII 



\H 2 (M) 

(iv) Next, let us deal with the last term 



V • q(T) - K V T Z 



wd z T 



V • q(T) - K V T ZZ dxdydz 



<C\\w\ 

Notice that from (61) 

lk|| L oc (n) < C 
By using (81) we obtain 



\8 Z T\ |V • q(T) - K V T ZZ \. 



o r 



AT(z,t) 



L°°(M) 



+ ||AT*| 



L°°(M) 



dz 



\\wh~ (a) <C[j h ||AT(^t)||V 2 2 (M) ||AT(^i)||^ 2 (M) dz + ||AT*|| L oo( M) 
By Cauchy-Schwarz inequality, (77) and Proposition 2 we get 

IIHk~ (n) < C (||AT*|| L oc (M) + UTII 1 / 2 |V • q(f) - K v f zz \^ 2 ) . 
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As a result of the above we reach 



/ 

Jn 

< C 



wd z T 



V • q(T) — K V T ZZ dxdydz 



(143) 



\T*\\ Hi{M) + \\T\\V 2 \V . q (T) - K V T ZZ \' 

x\d z T\ \V ■ q(T) - K V T ZZ \. 
Therefore, from the estimates (140)-(143) we have 
ld||T|| 2 



/2 



2 dt 



< C 



|V • q{T) - K V T Z , 



i + \Q\ 



V • q(T) - K V T Z 

+C\\f\\ 1 ' 2 \d z f\ \V-q(T)-K v T zz \i. 
By using Young's inequality we get 



H 4 (M) 



\d z T\ 



dt 

< C 



T\\ 2 + 



V • q(T) - K V T Z 



(144) 



i + iop + ini^^ + iT.i 4 



+ c 



1 + ||T|| 2 + |T Z 



Again, by Gronwall inequality we conclude 



where 



\\T(t)\\' 



K S (S,Q,T ,T* 



C 



V • q(T) - K v T zz \ 2 ds < K s (t,Q,T ,T*), 
1 + {K z (t, Q, To, T*)) 2 } t + K x {t, Q, To, T*) 



(145) 



(146) 



l^oll 2 + |Q| 2 + ||T1^ 4(M) +K z (t,Q,T ,T*) 



where Ki{t,Q,T ,T*) and K z (t,Q,T ,T*) arc as in (100) and (139), respectively. Since the strong 
solution is a weak solution, by Theorem 3, the strong solution is unique. □ 

5. Global Attractor 

In previous sections we have proved the existence and uniqueness of the weak and strong solution of 
the system (52)-(54). In this section we show the existence of the global attractor. Moreover, we give 
an upper bound, which are not necessarily optimal, for the dimension of the global attractor. Denote by 
T(t) = S(t)T the solution of the system (52)-(54) with initial data T . As a result of Theorem 3 and 
Theorem 5, one can show that 

f (t) = S(t)T e L 2 (fl) for all f £ T 2 (0), t > 0, 

and 

T(t) = S(t)T e V 2 for all T £ V 2 , t > 0. 
Since, in this section, we only consider the long time behavior of solutions of the system (52)-(54), by 
Theorem 3 and Theorem 5, we conclude that f(t) G £j£ c ((0, S], V 2 ) for every T(0) G L 2 {tt). As a result, 
one can easily show that 

d t TGT 2 oc ((0,S],V 2 '). (147) 

Theorem 6. Suppose that Q G L 2 (Q) and T* G H 4 (M). Then, there is a global attractor A C L 2 (VL) 
for the system (52)-(54). Moreover, A has finite HausdorfJ and fractal dimensions. 
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Proof. First, let us show that there are absorbing balls in L 2 (Q) and V2. Let T be the solution of the 



system (52)-(54) with initial datum T e L 2 (fi) and |T | = 
5, there is a to such that 

|T(i )| = T(i )+T* 



T + T> 



< 2p. 



< p. By Theorem 3 and Theorem 



From now on we assume that t > t . By taking the H' dual action to equation (53) with T, we obtain 

(d t f + V • q(f) - K v f zz ,f) + (v • VT + w9 z T + w • VT*, f ) = (<Q*,t) 
Applying (147) and Lions Lemma (cf. [31] Lemma 1.2. p. 260), we reach 



d t T,T) = - - 1 



Moreover, we have 



2 efe 

v- g (f)-^T zz ,r\ = ||f|| 2 . 



In addition, 



and 



v ■ VT + wd z T + v ■ VT*,T 



v ■ VT + w<9 z T + v • VT* 



Q* T dxdydz, 



T dxdydz; 



;q*,t; 

as long as the integrals on the right hand side make sense. Therefore, we have 
ld|T| 2 



2 dt 



+ \\T 



v ■ VT + wd z T + v ■ VT* 
By taking ip = T* in the weak formulation (70), we get 

[ f(t)T* dxdydz - [ f(t )T* dxdydz + [ a(f,T*) 
Jn Jn Jt 



T dxdydz = / Q* T dxdydz. 



+ 



t Jn 



v ■ VT + wd z T + v ■ VT* 



T* dxdydz = 



Q* T* dxdydz. 



t Jn 



It is equivalent to 
d 

dt.m 



4- f T(t)T* 
it J n W 



dxdydz + a(T,T*) + 



w-VT + w9 z T + v ■ VT* 



By adding (148) and (149), we obtain 

ld_ 
2~dl 



T* dxdydz = / Q* T* dccdydz. 



+ 



|T | 2 + 2 J f(t)T* dxdydz^j + ||f || 2 + o(f , T*) 
v VT + wd z f + v VT* (T + T*^J dxdydz 
= J Q* (T + T*) dxdydz. 



Notice that 



(148) 



(149) 
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(i) 
(ii) 



|T| Z + / 2TT* dxdydz = |T| 2 - / \T*\ 2 dxdydz; 
In Jn 



J 

Jn 



V-VT + wd z T + v • VT* 



(V + T*) dxdydz 



(iii) 



= / [v ■ VT + wd z T] T dxdydz; 
Jn 



T + T* ) dxdydz 



Q T dxdydz - a(T, T*) — \\T*\\ + a 
n J z=a 



TT* + — VT • VT* 

K v 



dxdy. 



Therefore, we get 



ld|T| 2 



2 dt 



+ ||T|| 2 + f [v ■ VT + wd z T] T dxdydz 
Jn 



J 

Jn 



Q T dxdydz + a 
n j z=o 



TT* + — VT • VT* 



dxdy. 



By integration by parts and (40), we obtain 

/ [v ■ VT + wd z T] dxdydz = 0. 
Jn 

By applying Cauchy-Schwarz inequality and (80), we reach 



/ 

Jz=0 



TT* + — VT • VT* 



dxdy 



< \\T*\\m(M) 



\T\\li(m) + -^||VT|| L 2 (M ) 



< Cl ,2 \\T 



l|r|| ff i ( n) + ^||VT|| H i (n) 



;C 2 1/2 (l + -|-)||T*|| ffl(M) ||T|: 



It is clear that 



/ 

Jn 



Q T dxdydz 



< \Q\ \T\. 



Therefore, by the above estimates (150)-(152), we obtain 



d\T[ 



C\\ 2 <C 2 a 2 [l + -^) \\T*\\ 2 H1{M) + 2\ 



dt \ K 

By Cauchy-Schwarz inequality and (74), we have 
d\T\ 2 1 



\T\. 



;/; +^\\Tr<C 2 a^l + f-) \\T*r H1{M) +2C \Q\ 2 , 
where Co is as in (74). Thus, again, by (74), we obtain 

« + ^|T| 2 <^(l + £) 112-11^+20,^ 



(150) 



(151) 



(152) 



(153) 
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By Gronwall Lemma, we get 



\T(t)\ 2 <\T(t )\ 2 e 2Co" + 2CoC 2 a^l + £) T \j, , , u , + 4^|0| 
As a result of the above, when t is large enough such that 



■»2|/-»|2 



■t 



|T(i )| 2 e 2C " < 2C C 2 a 2 1 



we have 



where 



In particular, 



\T(t)\ 2 <R a (T*,Q), 



njJ .Q)= WW (! + £;) II- 7 117'/', u, ^» Ql 



^2|/0|2 



(154) 
(155) 



/i 



limsup|T(^< 2C C 2 a 2 \l + j l|T*||^ 1(A/) +4C^|Q|^ 

In other words, when t is large enough we have 

\f(t)\ 2 = \T(t)-T*\ 2 <R a (T*,Q), 

where 



Ra(T*,Q) = 2R a (T*,Q) + 2\\T*\\ 2 L2 



(M)- 



(156) 



(157) 



where R a (T*, Q) is as in (155). Therefore, there is an absorbing ball in L 2 (f2) with radius R a (T*,Q) for 
system (52)-(54). 

Next, we show that there is an absorbing ball in V 2 . First, notice that from (153), we have 



J t t+r \\T( S )\\ 2 d S <2\T(t)\ + 



4Co C 2 a 2 [l + £) \\T*\\ 2 Hl{M) +8C 2 \Q\ 2 



Therefore, by (156), when t is large enough, we get 

rt+r 



J r \\T( S )\\ 2 d S <K r (r,Q,T*), 



where 



K r (r,Q,T*) = 2R a (T*,Q) 



4C CW (! + -£) \\T*\\ 2 HHM) + SC 2 \Q\ 2 



(158) 
(159) 



and R a (T*,Q) is as in (157). 

From the proof of Theorem 5 we recall the inequality (137): 



d_ 

dt 



T 2 dxdy + K h \VT z \ 2 + K V \T ZZ \ 2 + A|V • {H T VT Z )\ 2 + M |VT^| 2 + 



a 



z=0 



A|V • (H T VT)\ 2 + K h \VT\ 



dxdy 



<c(|Q| 2 + ||T*|| 2 , 1(M) + ||r|| 2 )+c(| 



t\\ 2 + \\t 



* 1 1 2 



\H 4 (M) 



) (\\T(z = 0)\\h (M) + \T z \ 2 ). 
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Thus, 



< c 



(' 



\ 2 + \\T*\\ 2 HHM) 



\T\\ 



C \\T\\ 2 + \\T 



*||2 



H4 ( M V \k~ v 



■H^ = 0)||| 2(M) + |T z 



By applying the uniform Gronwall inequality (cf. for example, [33], p. 89) and (158), we obtain, when t 
is large enough 

a 

where r > is fixed and 



|r,(t)| 2 + ^\\T(z - 0)||| 2(M) < R z (r,T*,Q), 



(160) 



R z (r,T*,Q) = C 



K r (r,T*,Q) 



+ \\T*\\ 2 HHM) r+\Q\'r + K r (r,T* 1 Q) 



(161) 



C 



xe 



K r (r,T*,Q) + \\T*\\ 2 Hi{M) r 



Let us recall the inequality (144): 
d 



V • q(T) - K V T Z , 



< C 



i + ior + ini^ (M) + m 



+ c 



i + ||ri| 2 + |T z 



Thus, 



dimi 



dt 



C 



1 + llTll 



ITII 



<c^i + |Q^ + ||r*|| 4 H4(M) + |r z 

Thanks the uniform Gronwall inequality (158) and (160), we obtain, when t is large enough 

\\f(t)\\<R v (r,T*,Q), 

where r > is fixed and 

-K r (r,Q,T*) 



(R v (r,T*,Q)Y = C 



r + \Q\ 2 r + \\T*\\% HM) r + R z {r, Q, T*) r 



(162) 



(163) 



C [r + R z {r,Q,T*f r + K r (r,Q,T*)] 



xe 

Therefore, we have shown that there is an absorbing ball B in V2 with radius R v (r, T*,Q). From the 
proofs of Theorems 3 and 5 we conclude that the operator S(t) is a compact operator. Following the 
standard procedure (cf., for example, [6], [7], [8], [14], [33] for details), one can prove that there is a 
global attractor 

A = n t>0 S(t)BcV 2 . 
Moreover, A is compact in L 2 (Q) due to the compact embedding of V2 in L 2 (fl). 

In addition to the compactness of the semi-group S(t) one can show its differentiability on A with 
respect to the initial data. Therefore, one can use the trace formula (cf. [6], [7], [33]) to get an upper 
bound for the dimension of the global attractor A. Let T be a given solution of the system (52)-(54) 
with T € A. Since it is on the global attractor A, T is a strong solution to the system (52)-(54). It is 
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clear that the first variation equations of the system (52)-(54) around T read: 
dtX = F'(f) X , 



0; d xX \ t= _ h = 0; §§ 



0; g( X ).n| r =0, 



(164) 
(165) 
(166) 



z=0 

X(x,y,z,0) = C, 

where x are the unknown perturbations about T with a given initial perturbation ( e L 2 (fl) . Moreover, 
here 

F'(T) X = - [v • q( X ) - K vXzz + u ■ V(T + T*) + v ■ V% - wd z f + wd zX 
where v = (vi,V2),w are as in (58)-(61), and 

e-Xx{x, y, £, t) + fx y {x, y, t) 



Ml 



-de- 



U2 



-I 



w 



e 2 + P 

° r £Xx(x, y, j, t) + fxyjx, y, g, t) 

-hJ-h t 2 + P 

z -fXx(x,y,£,t) +ex v (x,y,t,t) 

z 2 + P 

[' -fXx(x,y,^t) + e X y{x,y^,t) 

-hJ-h ' " c 2 + P 

'"" r e&x(x,V,Z,t) ~ foXx(x,y,£,t) 



d£dz, 



+ 



h J-h 
2/o/ 



e 2 + p 



d^dz, 

+ 



(167) 



(168) 



(169) 



{-fXx(x, y, £, i) + 6X5,(0;, 2/, £, *)) 



+ 



(e 2 + /2) 2 

" ' r eA x (x,y,^t) - foXx(x,y,t,t) 



d^d-q + 



z + h 



h J-h 



e 2 + P 



2/o/ 



7{-fXx{x,y,^,t) +ex v {x,y,^,t)) 



, d£dz 



(e 2 + /2) 2 

It is not difficult to show that the above, coupled second order elliptic and linear parabolic, system has 
a unique solution x(i). Moreover, for t > 0, 

X(t) G K 2 . 

For any positive integer m we consider the volume element \xi{t) A • • • A Xm(i)| A m L 2 (o) > we have the 
following trace formula (cf. [6], [7], [33]) 



^lxi(*) A---AXm(t)lAm^(n) 



Tr 



(P m {t) o F'(T(t)) o P m (t)) | Xl (t) A • • • A Xm(*)lA^(Q) . 



which gives 

\Xi(t) A • • • A Xm(i)lA-^(n) = |Ci A • • • A C m \\m L , {n) exp^ Tr (p ro (s) o F'(f (a)) o P m (s)) da, (170) 

where Xi( s ): " " ' 1 Xm(«) are the solutions of (164)-(169) corresponding to the initial data Ci, ■ ■ ■ , Cm, re- 
spectively. Tr (j 3 m (s) o F'(T(s)) o P m (s)j is the trace of the linear operator (j 3 m (s) o F'(T(s)) o P m (s)^ 

and P m (s) is the L 2 (£l) orthogonal projector onto the space spanned by {xi( s )i ' ' ' ,Xm(s)}- Thanks to 
(170) {xi( s ), ' ' ' ,Xm(s)} are linearly independent for every s > if and only if {£1, • • • , Cm} are lin- 
early independent. Hence, from now on we assume that {Ci,-- - , Cm} are linearly independent. Let 



30 



C. CAO, E.S. TITI, AND M. ZIANE 



{tpi(s), • • • , V'm(s)} be an L 2 (f2) orthonormal basis of the space spanned by {xi( s )> • ■ ■ j Xm{s)}- Notice 
that {V>i(s), • • • , tp m {s)} arc in V2 for s > 0. Thus we have 

m 

Tr (p m (s)oF'(T(s))oP m (s)) = £ (>(T (s))^(s), ^(s)) . 

3=1 

Observe that 

(f'(T ( S ))^-( S ),^(s)) = -||^(sf + ^ [«j • V(T + T*) + uJ^t] ^( s ) dzdydz, 

where, for j = 1,2, ••■ ,m, Uj(x,y,z,s) and wj(x,y,z,s) are as in (167)-(169), respectively, but with 
V^-(s) replacing %. Following the same steps that led the estimates (126), we have 



Uj --V(T + r*)^(a) 



and 



< C 



VT 



Wjd z T ipj(s) 
Recall that |^>j| = 1, for j = 1,2, • • ■ , m. Thus, 

itj • V(T + T*) + wd z T ipj(s) dxdydz 



<c(||t|| + ||ti* 1(j0 ) \i>m\ HMh 
\>l>A*)\ k WM\\ l - 



VT 



<C(||T|| + ||T*|| H1(M) ) \\il>i(s)\\+C 
By using Young's inequality and the above estimate, we have 

(F'(f(s))^(s), ^(s)) < -\UM\ 2 + C [ ||T( S )|| 2 + \\T*f Hl{M) + 

By (79), we have 

m 

ll^'( s ll 2 > Ai + • • • + A m > CX im 2 . 

3=1 

As a result, we obtain 

Tr (P m {s) o F'(f(s)) o P m ( s )) < C\\vr? + C 

Hence, 

-j\r(p m { S )oF'{f{s))oP m { s )} ds 

< -cx im 2 + c 

Therefore, by applying (84), we get 

lim sup sup sup 





6 




2" 


T 




VT 





\\T(a)r + \\T*\\ 2 H HM ) 





6 




2" 


T 




VT 









6 




2 


/' 


T 




VT 


ds 


Jo 











\\T{s)rds+\\T*r m{M) +- t 

- t J Tr (P m (s) o F'(T(*))) ° P ro (s) 



101 <i 



i = 1, 



, m 



< -CX im 2 + X 4 (T*,Q), 
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where 



K^T*, Q) = C(R a (Q,T*)f \l + \\T*\\ 2 H1(M) + |Q| 2 " . 



(171) 



In order to guarantee — C\\m 2 + K 4 (T*,Q) < we need to choose m large enough such that 



^ c ^ K 4 (T*,Q) y 



Therefore, the Hausdorff and fractal dimensions of the attractor A can be estimated by (cf. e.g., [8]) 
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